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Using the semiclassical Boltzmann transport theory, we analytically consider dc charge transport 
in gapless electron-hole (both chiral and non-chiral) systems in the presence of resistive scattering 
due to static disorder arising from random quenched impurities in the background. We obtain the 
dependence of the Boltzmann conductivity on carrier density and temperature for arbitrary band 
dispersion in arbitrary dimensionality assuming long-range (~ 1/r) Coulomb disorder and zero- 
range white noise disorder [~ 5(r)]. We establish that the temperature and the density dependence 
of the Boltzmann conductivity manifests scaling behaviors determining respectively the intrinsic 
semimetallic or the extrinsic metallic property of the gapless system. Our results apply equally 
well to both chiral and non-chiral gapless systems, and provide a qualitative understanding of the 
dependence of the Boltzmann conductivity on the band dispersion in arbitrary dimensionality. 


I. INTRODUCTION AND BACKGROUND 

Much interest has focused over the last 10 years on the 
electrical transport properties of electron (and hole) sys¬ 
tems where the band dispersion deviates from the stan¬ 
dard E(k) = h^k^/2m parabolic metallic band dispersion 
{E, k, m respectively being single-particle energy, wave 
vector, and effective mass). The most important exam¬ 
ple of such a system in nature, which has been studied 
extensively over the last 10 yearsii^, is graphene^ which 
has two-dimensional gapless linear relativistic Dirac-like 
chiral energy dispersion defined by E(k) = hvpk, with 
fc = |k| as the 2D wave number and vp as the graphene 
velocity defining its linear dispersion. 

Very recent interest in the literature^ has been on 3D 
Dirac-Weyl materials which have chiral relativistic linear 
gapless band dispersion in three dimensions (i.e., akin to 
‘3D graphene’). For our purpose the distinction between 
(degenerate) Dirac and (non-degenerate) Weyl systems 
is unimportant, and we will refer to the linearly dispers¬ 
ing gapless chiral systems as Dirac materials (with a cer¬ 
tain degeneracy factor ‘g’, which counts the number of 
carrier ‘species’ or ‘flavors’ in the system as reflected in 
the spin and/or valley degeneracy - e.g., for graphene 
g = i and for some of the proposed 3D Dirac materials 
g = 24), completely ignoring any inter-species scattering 
(e.g., scattering between valleys in the band structure) 
and assuming each carrier species to be contributing in¬ 
dependently (and equally) to charge transport. (The 
neglect of inter-species (often, inter-valley) scattering is 
an excellent approximation in most experimental situa¬ 
tions since such scattering involves very large momentum 
transfer.) A number of recent theoretical publications 
in the literature study transport properties in 3D Dirac 
systems^^— . The number of publications on 2D graphene 
is huge (> 10,000) and we can only cite some relevant 
review articlesii^. 

Our current work presented in this paper asks the fol¬ 


lowing question: What can we say theoretically about 
transport properties of a gapless electron-hole system 
with arbitrary energy band dispersion (and in arbitrary 
dimensions)? Of course, the general theory for such a 
completely general problem is intractable, and we fo¬ 
cus on the semiclassical transport properties within the 
Boltzmann theory using the relaxation time approxima¬ 
tion. We also restrict ourselves to the dc electrical con¬ 
ductivity limited only by disorder scattering where the 
disorder arises from random quenched impurities in the 
material with electron-impurity scattering potential be¬ 
ing either the long-ranged Coulomb potential (assuming 
to be the usual 1 /r-type in all dimensions where r is the 
electron-impurity spatial separation) or the zero-range 
(5-function potential. We consider both screened and un¬ 
screened limits of the Coulomb disorder, and obtain com¬ 
pletely analytical results for the dc charge conductivity in 
all dimensions with our results agreeing with the known 
results in two and three dimensional gapless Dirac sys¬ 
tems. It turns out that our analytical results obtained in 
this work are very general, and apply to both chiral and 
non-chiral electron-hole gapless systems although most 
existing gapless systems (e.g. graphene) happen to be 
chiral in nature, obeying some effective form of massless 
Dirac-Weyl equations at long wavelength. 

Our motivation for our theoretical work is obvious. 
Civen the great current interest in gapless linearly dis¬ 
persing Dirac systems, it makes sense to obtain results 
for both linear and quadratic energy dispersions in ar¬ 
bitrary dimensions in order to develop intuition about 
the dispersion-dependence of transport properties. Note 
that we restrict the system dimensionality to be greater 
than 1 (i.e., d > 1, where d is the system dimensionality) 
throughout this paper. If theoretical results can actu¬ 
ally be obtained in arbitrary dimensions ((?>!) for ar¬ 
bitrary band dispersions, obviously that would provide a 
more complete picture for the qualitative nature of trans¬ 
port properties in electronic materials, and doing this is 
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therefore our primary motivation. But there is actually 
a second subtle (indirect) reason arising from quantum 
criticality for our theoretical work, which we discuss be¬ 
low. 

It has been known for almost 30 years, since the semi¬ 
nal work of Fradkinii, that a system of gapless electrons, 
where conduction (electron) and valence (hole) bands 
touch at a finite number of isolated points (so-called 
Dirac points where the system is a perfect semimetal 
with no free electrons or holes at zero temperature) in 
momentum space manifesting linear band dispersion, has 
a disorder-induced T = 0 quantum phase transition for 
d > 2 (with dc = 2 being the critical dimension). The 
transition manifests itself as the density of states (DOS), 
D{E), becoming finite at FI = 0 (taking the Dirac point 
to be at = 0 with no loss of generality as our notation 
for zero energy) for infinitesimal disorder strength (finite 
critical disorder strength) for d = 2 (d > 2). For lin¬ 
early dispersing carriers, of course, D{E) goes as E {E"^) 
in d = 2 (d = 3) in a clean noninteracting system with 
D{E = 0) = 0 at the Dirac point in both cases. The DOS 
serves as the nominal order parameter for this quantum 
phase transition, being zero (finite) for disorder strength 
below (above) the critical point (with the critical disor¬ 
der being zero for d = 2). Although the original result 
was obtained within an 1/A^-expansion with N = g be¬ 
ing the number of fermion species (or the electron degen¬ 
eracy factor), this disorder-driven quantum phase tran¬ 
sition has now been verified using many disparate the¬ 
oretical techniques including one-loop^^ and two-loopi^ 
renormalization group (RG) calculations, exact numeri¬ 
cal simulations^ d^d^ ^ self-consistent Born approximation 
calculations^i^, and general scaling theoriesi^. In fact, 
very recent workii indicates that this disorder-induced 
quantum phase transition even survives some finite in¬ 
teraction effects leading to rich multicritical behavior. 

We emphasize five important caveats with respect to 
this quantum phase transition: (i) it occurs in the un¬ 
doped intrinsic situation when the chemical potential is 
precisely at the Dirac point {E = 0) with no free elec¬ 
trons/holes present in the system (indeed for any finite 
doping, the chemical potential would shift from the Dirac 
point into the conduction or the valence band leading to 
D{E = Ep) at the fermi level to be finite any way even 
in the clean system indicating the system to be a ’metal’ 
with finite DOS); (ii) the transition has only been theo¬ 
retically established for short-ranged (e.g., Gaussian) dis¬ 
order; (iii) the so-called Griffiths physics associated with 
rare disorder regions may affect the transition leading to 
band tailing corrections even at E = which could lead 
to finite DOS for all disorder in all dimensions (thus ren¬ 
dering the transition into a crossover); (iv) the physics of 
disorder-induced electron-hole puddle formation^ which 
is known to be of great importanc o^^i^° and relevance^i^ 
in graphene, is ignored (indeed disorder induced puddles 
produce to macroscopic density inhomogeneity in the sys¬ 
tem with considerable spatial fluctuations in the chemi¬ 
cal potential which lead to the local DOS being always 


a fluctuating finite quantity throughout the system for 
any finite disorder even when the average Fermi level is 
at E = 0] (v) all quantum localization effects are ignored 
consistent with the assumption of no intervally scattering 
in the system^^. 

The disorder-induced ‘semimetal-to-metal’ transition 
(with the DOS increasing from zero in the semimetal to 
finite in the metal at the Dirac point) discussed above 
provides an important (albeit indirect) context for our 
current work. The Boltzmann conductivity at T = 0 for 
an electron system is given by cr oc D{Ep)t[Ep), and 
obviously the DOS plays a decisive role in directly de¬ 
termining the conductivity through the D{Ep) factor as 
well as indirectly through the relaxation time r which 
itself, of course, depends on the DOS (and also on the 
disorder potential). But, the DOS, D{E), for an ar¬ 
bitrary band dispersion defined by E(k) ^ goes as 
D ~ E‘^/°‘-^ in d-dimensions (e.g., for linear dispersion, 
a = 1, D E, E^ in d = 2, 3, and for parabolic disper¬ 
sion, a = 2, D ^ E°, E^/^ in d = 2, 3, respectively, as 
is well-known), and thus specific combinations of {a,d) 
provide finite or vanishing DOS at E = 0. Thus, for 
d = 2 = a (i.e., the ordinary parabolic metallic elec¬ 
tron system), the DOS is an energy-independent constant 
which is finite at FI = 0 whereas for d = 3 = a the same 
is true in 3D systems. In fact, for any d = a {d > a) the 
DOS is finite (zero) at FI = 0 in arbitrary dimensions for 
the clean noninteracting system, thus making a = 2 (3) 
special for d = 2 (3) respectively in the consideration of 
the putative semimetal-metal theoretical consideration. 
We note that the critical dimensionality for the disorder- 
driven (through the DOS) semimetal-metal transition is 
precisely dc = 2a^ showing the coupled role of d and 
a in this context. Thus, we can compare the effect of 
having finite versus vanishing DOS on transport prop¬ 
erties by tuning the band dispersion (i.e., a) and the 
dimensionality (i.e., d). This should provide some quali¬ 
tative insight into the behavior of the charge conductivity 
through the semimetal (“vanishing DOS”) to metal (“fi¬ 
nite DOS”), albeit by tuning a (rather than by disorder). 
We are also in a position to make statements about how 
the nature of disorder (long versus short range) affects 
the conductivity in the semimetallic versus the metallic 
phase. 

Thus, the theoretical issue we address in the current 
paper is complementary to that considered by Fradkin 
and others in the context of disorder-driven quantum 
criticality2rJ^ - we ask how the Boltzmann conductiv¬ 
ity behaves as a function of the band dispersion param¬ 
eter (a) in different dimensions (d > 1) for long- and 
short-range disorder assuming that the DOS of the gap¬ 
less system is determined only by band dispersion (and 
dimensionality), and not by disorder (with disorder enter¬ 
ing the calculated Boltzmann conductivity only through 
the transport relaxation time t). Since the clean sys¬ 
tem DOS depends on a and d (and may be finite or zero 
at the Dirac point, depending on the precise values of 
a and d), we anticipate deriving general relations be- 
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tween the conductivity of the gapless system (within the 
Boltzmann theory where the DOS is not renormalized by 
disorder) and values of a and d (which in turn determine 
whether the clean system is a semimetal with vanishing 
DOS or a metal with finite DOS). We are therefore estab¬ 
lishing a connection between the clean system DOS and 
the semiclassical conductivity, and not exploring quan¬ 
tum critical physics. From the perspective of quantum 
criticality, the undoped gapless system has a critical di¬ 
mensionality of dc = 2a, implying that for a = 1 (2), 
dc = 2 (4) so that for linear (quadratic) dispersion, the 
system can only have a disorder-induced semimetal-to- 
metal transition with the DOS as the order parameter 
for three (five) dimensions (or above). For a = 3, dc = 6, 
and thus the disorder-driven semimetal-to-metal transi¬ 
tion considered originally by Fradkinii becomes trivial 
in this case. Clearly, graphene (d = 2 and a = 1) is 
the marginal case for linearly dispersing systems whereas 
d = 4 is the marginal case for parabolic dispersion case 
(a = 2). This makes a < 3/2 (1) be the interesting 
nontrivial band dispersion case for 3D (2D) gapless sys¬ 
tems. We refer to clean systems with vanishing (finite) 
density of states at the Dirac (i.e. band touching) point 
{E = 0) as semimetal (metal) throughout this paper, and 
theoretically investigate the resultant Boltzmann conduc¬ 
tivity arising from short-range and long-range impurity 
scattering. We are thus asking a question different from 
(and complementary to) that discussed in Refs. [Ml- 

We note that the order parameter controlling the ac¬ 
tual disorder-driven quantum phase transitionii^— is 
the DOS, being zero (finite) in the semimetal (metal) 
phase, and therefore any connection between this quan¬ 
tum phase transition and the corresponding behavior of 
the charge conductivity, i.e., the semimetal (metal) be¬ 
ing defined by vanishing (finite) conductivity, can only be 
established numerically. In particular, within the Boltz¬ 
mann theory, D{Ep = 0) = 0 does not necessarily imply 
a vanishing conductivity since the transport relaxation 
time could diverge a.t Ep = 0, and conversely, one could 
have, at least as a matter of principle, a vanishing con¬ 
ductivity for finite values of D{Ep) since the relaxation 
time could vanish! One of our goals in the current work is 
precisely understanding this dichotomy within the Boltz¬ 
mann theory by investigating various combinations of d 
and a to see whether it is indeed always true that the 
DOS and conductivity manifest monolithic behavior with 
the vanishing (finite) values of D{Ep = 0) being neces¬ 
sary and sufficient for vanishing (finite) conductivity in 
all situations. We find the answer to this question within 
the Boltzmann theory to be a definitive negative. 

The question of how the Dirac point charge conductiv¬ 
ity behaves through the semimetal-metal quantum phase 
transition has been somewhat controversial^^— In 
particular, although the conductivity on the metallic side 
of the transition is always claimed to be finite by virtue of 
the hnite DOS, the conductivity on the semimetal side 
(with vanishing DOS) has been found to be vanishing 
(finite) for short-range (long-range) disorder within the 


self-consistent-Born-approximation (SCBA) for 3D Dirac 
systems (d = 3 and a = 1) in ref. Q. The direct numeri¬ 
cal calculation^ finds finite conductance at A = 0 in both 
phases, but infers that the conductivity actually vanishes 
(remains finite) in the semimetal (metal). Our earlier 
Boltzmann theory (using the standard Born approxima¬ 
tion for obtaining the relaxation time) gives vanishing (fi¬ 
nite) conductivity at the Dirac point for both 2D2i and 
3D^ Dirac systems using long-range (short-range) disor¬ 
der in contrast to the SCBA results. Thus, the physics 
for long-range disorder, which is of considerable exper¬ 
imental relevance because of the invariable presence of 
random charged impurities in electronic materials, re¬ 
mains open as does the question of whether the inter¬ 
esting semimetal-to-metal transition (defined as a tran¬ 
sition from vanishing to finite DOS) at the Dirac point 
in clean 3D Dirac materials can be distinguished exper¬ 
imentally (at least as a matter of principle) by studying 
the electrical conductivity through the transition. 

The rest of the paper is organized as follows. In sec¬ 
tion II, which gives the main content of our work, we 
present our Boltzmann transport theory providing ana¬ 
lytical results for the conductivity for various values of 
a and d. We provide a discussion of our results along 
with a conclusion in section III. An appendix discussing 
the temperature dependence of screening completes the 
presentation. 


II. THEORY AND RESULTS 

Although our primary interest is understanding the be¬ 
havior of the T = 0 Boltzmann conductivity at the Dirac 
point for the undoped Dirac material (i.e., the chemi¬ 
cal potential at the Dirac point, kp = Ep = 0, with no 
free electrons or holes in the system) as a function of di¬ 
mensionality (d) and band dispersion (a) for long-range 
and short-range disorder, we will consider, for the sake of 
generality, a system with a finite doping {kp,Ep ^ 0) so 
that our results have general applicability to experimen¬ 
tal systems which can never be ideally intrinsic with the 
chemical potential precisely at the Dirac point. It is easy 
to figure out the intrinsic Dirac point behavior by taking 
the appropriate kp,Ep —>■ 0 limit of our general density- 
dependent results. We present our analytical results in 
three different subsections with IIA, B, C giving respec¬ 
tively results for arbitrary dimensions (d > I), d = 3 
(and 2), and finally (IIC) the results for hnite tempera¬ 
tures (with section IIA and B being restricted to T = 0). 
We use the word “Dirac point” to signify the E = 0 band 
touching point for all values of a (i.e., arbitrary energy 
dispersion) although the word Dirac point strictly applies 
only to the linear chiral band dispersion case (a = 1). As 
mentioned already, our results remain valid for all gapless 
systems independent of their chirality. 
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A. Boltzmann Conductivity with E = Ak°‘ in 
arbitrary dimension (d > 1) 

We first consider Dirac materials in arbitrary dimen¬ 
sions {d > 1). For a clean system with an energy disper¬ 
sion 


Ek = A\kr, ( 1 ) 


(1) Screened Coulomb potential - within the Thomas- 
Fermi (TF) approximation the screened Coulomb poten¬ 
tial can be expressed as 


Vd{q) 


+ qs~^) 


( 7 ) 


and the screening wave vector in d-dimension is given 
by (gs = qTF in the TF-approximation) 


the density of states in a d-dimensional space is given by 

where r(a;) is the gamma function, and the carrier den¬ 
sity (n), which is defined as particle number per unit 
d-dimensional volume (with only d = 2, 3 being of phys¬ 
ical interest), is related to the Fermi wave vector {kp) 
as 


"" 2‘^^^/2r(d/2 + l)^^’ 

where g is the total degeneracy (i.e., number of fermion 
flavors including both valley and spin). Thus, we have 
kp oc 

The Boltzmann charge conductivity can be expressed 
as 

. = fj/ dED^(EKr,, (-§fe) . (4) 

where Vk = dEk/dk is the velocity of the carrier, and fp^ 
is the Fermi distribution function. The scattering time 
Tp^, in d-dimension (ignoring intervally scattering in the 
system) is given by 

I _ Hi Ef! 

X Jd‘^-^n\v{q)\'^Un){i-cos9), (5) 

where Ui is the d-dimensional impurity density, Ep = 
Akp is the Fermi energy, d‘^“^D is the (d—1) dimensional 
solid angle element, /a(f2) is the wavefunction overlap 
factor arising from the chirality of a system, 0 is the 
scattering angle between k and k', and g = 2fcsin(d/2). 
Note that for a nonchiral system /a(D) = 1 in Eq. ([5]). 

The Coulomb interaction in the wave vector space is 
given by the appropriate d-dimensional Fourier transform 
of the Coulomb interaction V{r) = i.e., 


Vd{q) 






d-l 


( 6 ) 


where k is the background lattice dielectric constant of 
the material. Note that we take the Coulomb potential 
to be 1/r in all dimensions (d > 1) where ‘r’ is the d- 
dimensional spatial separation. 


qt^ = qT~P = dyep). ( 8 ) 

From Eqs. @, ®, and ® we can calculate the density 
dependent scattering time and conductivity at zero tem¬ 
perature. In the strong screening limit (gg ^ kp) the 
scattering time becomes 

T oc Ep°‘~^ oc oc (9) 

and the Boltzmann conductivity becomes 

CT oc oc oc (10) 

(2) Unscreened Coulomb potential or weak screening 
limit {kp ^ qs )— in this case the scattering time is cal¬ 
culated to be 

r oc o, fcd- 2 +a ^ „l+(a- 2 )/d (^^ 

and the conductivity to be 

a oc 42d-4)/a+2 ^ fcM-4+2„ ^ ^2+(2a-4)/d ( 13 ) 

(3) Density independent screening wave vector, i.e. 
qs = constant - in the strong screening limit g^ 3> kp, 
we have 

T oc oc kp~‘^ oc 

a oc 4 ”^''“ oc 4-2 „ 2 (a-l)/d^ (i 3 ( 

and in the weak screening limit {kp ^ g^) the results are 
the same as the results for unscreened Coulomb potential. 

(4) Short range potential, i.e., Vq = constant - in this 
case the results are the same as Eq. dni) 

r oc 4 ”^^°^ oc kp~'^ oc 

a oc oc 4“^ oc (14) 

Thus, for n —> 0 the conductivity is only finite for the 
short range potential and a = 1 regardless of the system 
dimensionality. (Note that here d > 1 throughout.) For 
all long-ranged disorder, screened or not, the semimetal 
phase (n = 0) has vanishing Boltzmann conductivity for 
all disorder (except for the density-independent screen¬ 
ing, which is exactly equivalent to short-range disorder 
since kp = 0 aX the Dirac point). 

In the next subsections, we use the general d- 
dimensional analytical results derived above to obtain 
specific Boltzmann conductivity formula in physical di¬ 
mensions (d = 2, 3) for arbitrary band dispersions and 
for different types (e.g. the long- and short-range disor¬ 
der potential) of resistive impurity scattering in gapless 
electron-hole systems. 
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B. Conductivity in 3D (2D) with arbitrary band 
dispersion 


In this subsection we obtain the density dependence of 
the Boltzmann conductivity for the experimentally rele¬ 
vant 3D (and 2D) gapless systems. For screened Coulomb 
disorder we have {d = 3) 


(D(k,k')) 


47re^ 1 

K + g2 ’ 


(15) 


where Qs is the screening wave vector in 3D. For = 0, 
we get the usual unscreened Coulomb potential in the 
wave vector space going as 1/g^ indicating the long-range 
l/r behavior of the Coulomb interaction. From Eq. (21) 
we have the density of states per electron flavor or species, 
Dct{E), in 3D as (with E ^ k°‘ for the gapless electron- 
hole band dispersion) 


and 

fc2a+2 „2(a+l)/3 

<7 OC — - OC - 

Hi 


( 22 ) 


Thus, for unscreened disorder the conductivity is a non 
universal function of density (i.e., dependent on a), but 
it nevertheless vanishes in the n —>■ 0 semimetal limit for 
all values of a (and arbitrary chirality). 

For short range disorder V{q) is a constant and we 
have 


T OC — - OC -, 

m Hi 

giving 

a OC — -OC- 

m rii 


(23) 


(24) 


Do^iE) = 


1 

27r2 aA^/'^ 


(16) 


and within the 3D TF approximation we have for the 
screening wave vector from Eq. ([5]) 

9s = 9tf =- —Da{Ep) OC Ep . (17) 


In this case the conductivity becomes finite for only a = 1 
as n —>■ 0, and the semimetal Dirac point manifests a 
finite conductivity in the Boltzmann theory for all values 
of the disorder strength (and chirality). 

For a 2D system we can calculate the density depen¬ 
dent conductivity also, obtaining 

(1) strong screening Coulomb disorder [q^ ^ kp) 


Since Ep oc we have q^ oc where kp is the 

Fermi wave vector which is related to the carrier density 
through n = gkp/()Tr^. 

The Boltzmann conductivity at zero temperature is 
given by 


a = e^^gDa[Ep)T (18) 

where vp = dE^/dk oc kp~^. In the strong screening 
limit, qs ^ kp, we have 


2 1.3-“ 

Qs 

T CX - OC - OC 

rii m 


n. 


(19) 


Combining with Eq. (1181) we have for three dimensions 


n 

a oc — oc 


4/3 


m 


rii 


( 20 ) 


for all a. Thus the conductivity vanishes at the Dirac 
point (n —> 0). This result is universal and is independent 
of the energy dispersion (i.e., a) and chirality (i.e. the 
functional form of / in Eq. ([5|)). Thus, the Boltzmann 
theory already predicts zero conductivity at the Dirac 
point for screened Coulomb disorder, and the semimetal 
state is preserved independent of disorder strength (i.e., 
value of rii). This result applies for arbitrary band dis¬ 
persion and chirality in gapless electronic systems at the 
Dirac point. 

For unscreened Coulomb potential or very weak 
screened potential (kp ^ qs) we have 


kjp 

a oc — oc n. 


(25) 


(2) unscreened or weak screening Coulomb disorder 
(qs < kp) 


cr CX —^ CX — 

rii n. 


(26) 


(3) short-range disorder 


a CX 


oc 


m 


rii 


(27) 


Again, for linear Dirac dispersion (a = 1), the 2D Dirac 
point conductivity is finite for the semimetal for short- 
range disorder only, vanishing for screened or unscreened 
disorder. These 2D results also apply for all gapless sys¬ 
tems independent of their chirality. 


C. Temperature dependent conductivity arising 
from energy averaging in the semimetal 

Before considering the temperature dependent conduc¬ 
tivity we first calculate the activated carrier density at 
finite temperatures. For a gapless intrinsic d-dimensional 
system, in which conduction and valence bands meet at 
the Dirac point and therefore the chemical potential is lo¬ 
cated at the Dirac point (/r = 0), the activated electron 
density is given by 


T (X — - oc - 

rii rii 


( 21 ) 


n(T) = gj 


d‘^k 1 

(2Tr)‘^ -|- 1 ’ 


(28) 
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where g is the total degeneracy and /3 = l/ksT. We can 
rewrite this equation as 

3 1 d 

’ 2'^7r'^/2 T{d/2 + 1) Ad/^a Jo + 1 

= C{kBTf^°‘, (29) 

where C is a constant. Thus the activated electron 
density is proportional to In 2D (3D) the ther¬ 

mally activated density increases linearly with tempera¬ 
ture for quadratic (cubic) dispersion E = {E = Ak^) 
whereas for linear band dispersion, the activated density 
increases as (T^) for 2D (3D) systems. 

The conductivity Eq. can be written as 

a = B dEE^+^r{E) (30) 

where S is a constant. Assuming that t oc E^ and is 
independent of the temperature, the energy averaging 
of the scattering time gives the following temperature 
dependent Boltzmann conductivity in the n —> 0 limit 
(i.e., intrinsic undoped systems) 

ao^T^+^+^. (31) 

The energy dependent scattering time is given in Eqs. 
(HID, (USD, and da for different situations, and substi¬ 
tuting the corresponding /3 in Eq. dSID one obtains the 
temperature dependent Boltzmann conductivity for the 
semimetal (i.e., n = 0 ) for different types of disorder. 

Writing the temperature dependent Boltzmann con¬ 
ductivity in the semimetal (i.e. Ep = 0) as 

a{T) - T^, (32) 

we get 

d-2 ^ , 

7 —-h 1 -I- /3, (33) 

a 

where (3 is the exponent defining the energy dependence 
of the transport relaxation time (r ^ E^). We note that 
although Eqs. (15^ and (1331) give the main temperature 
dependence of the Boltzmann conductivity, we have ig¬ 
nored the (presumably weak) T-dependence arising from 
screening (see Appendix A) as we assume qs to be tem¬ 
perature independent in deriving Eq. (IMp . 

Eor specificity, we note that in 3D linear dispersion 
Dirac semimetal (d = 3 and a = 1) we get 7 = 4, 4, 0 
for screened Coulomb, unscreened Coulomb, and short- 
range disorder, respectively (the corresponding 7 = 2, 3, 
1 for 3D chiral quadratic dispersion, d = 3 and a = 2). 
Thus, the temperature dependent semimetallic {Ep = 0) 
conductivity of the gapless system, in principle, provides 
direct information about both the energy dispersion (and 
hence the noninteracting DOS) and the disorder poten¬ 
tial (i.e., long-range or short-range). We note here also 
that the scaling of the Boltzmann conductivity, a{Ep) 
at T = 0, on the Eermi energy Ep, given by tr ^ Ep 


with the exponent S (as given above in Eqs. ® - (fT4l) h 
follows exactly the same scaling law as thermal scaling, 
i.e., d = 7 as will be discussed further later in section 
III. (The equality d = 7 holds quite generally up to loga¬ 
rithmic accuracy provided that screening itself is weakly 
temperature-dependent and the main temperature de¬ 
pendence arises from the thermal averaging.) 

In particular, let us consider d = a = 3 whence the 3D 
DOS is finite in the semimetal phase (instead of vanish¬ 
ing at the Dirac point, EiT- = 0, as in the linear dispersion 
situation). Eor d = a = 3, we have 7 = 4/3, 8/3, 4/3 
for screened Coulomb, unscreened Coulomb, and short- 
range disorder, respectively. This implies that, although 
the DOS itself is finite for d = a = 3 case, the correspond¬ 
ing Dirac point Boltzmann conductivity always vanishes 
for all disorder types keeping the system a semimetal (as 
far as the conductivity goes) independent of the fact the 
DOS itself is finite! By contrast, d = 3 and a = 1 case has 
7 = 4 (screened Coulomb), 4 (unscreened Coulomb), 0 
(short-ranged), implying that the Dirac point Boltzmann 
conductivity vanishes for Coulomb disorder (screened or 
unscreened), but is finite for short-ranged disorder al¬ 
though the DOS at E = 0 vanishes in this case. Thus, 
we actually get the intuitively unexpected result that in 
3D clean systems, in the presence of short-range disorder, 
the ‘semimetallic’ state with vanishing DOS at E = 0 
(i.e. the a = 1, d = 3 case) has finite conductivity 
whereas the ‘metallic’ state with finite DOS at E = 0 
(i.e., a = 3, d = 3 case) has vanishing conductivity! 
Thus, at least for the clean system, where the DOS can 
be tuned by the dispersion exponent a, there is no con¬ 
nection between the DOS and the Boltzmann conductiv¬ 
ity. Given that the d = 3, a = 1 case has vanishing DOS 
at Ep = 0 whereas the d = 3, a = 3 case has finite DOS 
at Ep = 0, we conclude that the temperature dependent 
intrinsic semimetallic conductivity (i.e., the actual value 
of the exponent 7 ) could, in principle, distinguish the sit¬ 
uation of finite DOS from the vanishing DOS in a gapless 
system, but not whether the conductivity itself is finite 
or not. All the temperature-dependent results derived 
in this sub-section apply both to chiral and non-chiral 
gapless systems. 


III. DISCUSSION AND CONCLUSION 

Our main theoretical finding is that the semiclas- 
sical Boltzmann conductivity always vanishes at the 
semimetallic “Dirac point” (i.e., Ep = 0) for all clean 
gapless systems in all dimensionalities for all Coulomb 
disorder (screened or unscreened) and short-ranged dis¬ 
order except for the short-range uncorrelated white-noise 
disorder in the linear dispersion case (a = 1 ) which pro¬ 
vides finite dc conductivity in all dimensions for all disor¬ 
der strength. We have obtained general analytical results 
for the conductivity, a{Ep,T,a,d]ni), as a function of 
Fermi energy (or equivalently carrier density), temper¬ 
ature, band dispersion {E ^ A:“), system dimensional- 
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ity, and disorder strength or type with Ep = 0 denoting 
the “semimetallic” intrinsic phase and Ep ^ 0 being the 
generic “metallic” extrinsic phase with free carriers. In 
the presence of (screened or unscreened) Coulomb dis¬ 
order and short-ranged disorder with a ^ 1, we find 
cr(T = 0; Ep = 0) = 0 for all values of d as long as 
rii ^ 0 whereas for short-range disorder in the linear dis¬ 
persion situation we find a{T = 0; Ep = 0) 0 for all 

d and a = 1. Our qualitative results with respect to 
conductivity (i.e., whether a is zero or finite) thus can¬ 
not distinguish between zero of finite DOS phases in the 
3D Dirac system (with a = 1) dX Ep = 0 since both 
situations (vanishing or finite DOS at E = Ep = 0) will 
give either zero conductivity (Coulomb disorder) or finite 
conductivity (short-range disorder) at T = 0. Thus, the 
Boltzmann transport theory does not directly manifest 
any direct or indirect effect of the semimetal-to-metal 
quantum phase transition at the Dirac point. It is, how¬ 
ever, worth while to mention that the pure Dirac linear 
dispersion case (a = 1 ) is special (at least with respect to 
short-range disorder) since it is the only situation which 
allows for a finite dc conductivity (in both d = 2 and 
3) at Ep = 0 in spite of the DOS being zero. If we use 
Boltzmann conductivity for Ep = 0 as the ‘order param¬ 
eter’ (rather than the DOS itself), we conclude that the 
disorder-driven critical transition discussed in Refs. 0- 
[Tsjl should not be called a semimetal-metal transition 
simply because the DOS changes from being zero below 
to being finite above the transition since the expression 
‘metal’ typically implies a finite T = 0 conductivity, and 
there is no evidence whatsoever that this is true in the 
putative ‘metallic’ phase for long-range disorder as T and 
Ep go to zero. A recent direct numerical worb^ investi¬ 
gating the conductance for disordered 3D Dirac systems 
explicitly considered short-range disorder. 

Our detailed analytical scaling results for a{Ep,T), 
however, suggest some qualitative possibilities for dis¬ 
tinguishing between semimetallic {D = 0 at Ep) and 
metallic (D ^ 0 a.t Ep) phases through the power law 
dependence of the conductivity on Ep (or equivalently 
density n) and/or temperature T. We find (see section 
II for details) 


a{T,Ep = 0) ~ T'^, 

(34) 

a{Ep,T = 0)^E%, 

(35) 


where the scaling exponents 7 and 6 depend on a, d, 
and disorder. While 7 is defined in Eqs. (15^ and (1551) . 
we provide S below for various situations (by collecting 
together the results from section II) 

{ 2{d —l)/a screened Coulomb disorder 

-I- 2 unscreened Coulomb disorder . 

a. 

2 — 2 /a, short-range disorder 

(36) 

We note that for density-independent screening (i.e., as¬ 
suming that the screening wave vector qg is independent 



FIG. 1. The figure shows the metal-semimetal generic 
phase diagram defined by the Boltzmann condnctivity with 
the high-temperature (T > ITfI) regime manifesting the 
semimetallic scaling cr ~ T'*' whereas the low-temperature 
(T < |Tf|) regime manifesting the metallic scaling cr ~ Ep 
(with 7 = <5 in the simple theory). All conductivities except 
for linear dispersion with short-range disorder gives vanishing 
conductivity at the Dirac point at T = 0. The explicit ana¬ 
lytical values of the scaling exponents 7 = <5 are given in the 
text. 

of kp and is a constant), the strongly screened Coulomb 
disorder gives the same results as the short-ranged dis¬ 
order. Comparing Eqs. and (1551) . and putting in 
appropriate /3 values, we conclude (5 = 7 , which is an 
important verification of our theory since we expect this 
exponent identity based on general scaling arguments. 

We note that the scaling exponents 7 = i5 in the Boltz¬ 
mann theory do in general depend on both a and d (as 
well as on /3, with r ~ E^, which is determined by the 
type of disorder) except for the pure short-ranged dis¬ 
order which is independent of dimensionality. Thus, we 
have the following identity valid for all disorder (provided 
the temperature dependence of screening is negligible and 
the temperature dependence in the conductivity is arising 
entirely from energy averaging in the Boltzmann theory) 

7 = A (37) 

The “semimetallic” behavior in conductivity manifests 
itself ior T > Tp = Ep/kp with a ~ T'*' and the “metal¬ 
lic” behavior manifests itself for T < Tp with cr ~ Ep, 
and (5 = 7 . This physics is schematically depicted in 
Fig. 1 which also serves to summarize our theoretical re¬ 
sults showing that the “high-temperature” (i.e. T > Tp 
which is also the low-density) regime corresponds to the 
semimetal and the “low-temperature” (i.e., T < Tp, 
low density) regime corresponds to the “metallic” regime 
with a finite T = 0 value of conductivity. 

We note that the Boltzmann theory predicts cr(T = 
Q,Ep = 0) = 0 even for a = 3 (2) in d = 3 (2) dimen¬ 
sions for all disorder models (including short-range dis¬ 
order) although the corresponding 3D (for E ^ k^) and 
2D (for E ^ k"^) DOS is finite at Ep = 0. Thus, within 
the Boltzmann transport theory at least, having finite 
DOS at Ep = 0 is not synonymous with a finite metallic 





conductivity at Ep = 0, and one should be careful in 
automatically identifying metals (with finite conductiv¬ 
ity) and semimetals (with zero conductivity) with having 
finite or zero DOS respectively! 

We note that to the extent exact numerical work^ 
and SCBAi find the disorder-driven semimetal-to-metal 
quantum phase transition to be directly connected with 
the Dirac point conductivity increasing from zero to a fi¬ 
nite value at a finite disorder (also as shown in Ref. ini), 
the disorder model used in these theories is short-ranged 
disorder which provides a nongeneric answer even in the 
Boltzmann theory of ours where in d = 2 and 3 we find 
the linearly dispersing Dirac system to always have a fi¬ 
nite conductivity even in the T = 0, Ep = 0 limit in the 
presence of only short-range disorder. For long-range dis¬ 
order, the conductivity vanishes for Ep = 0, and indeed 
the SCBA theory also no longer gives results consistent 
with the putative quantum phase transitioni^ The issue 
of quantum criticality associated with long-range disor¬ 
der remains an open question at this stage. 

Our work indicates that the semimetal-metal transi¬ 
tion, even if it is theoretically allowed, will be very dif¬ 
ficult, if not impossible, to observe experimentally. Any 
real 3D linearly dispersing (a = 1) Dirac system will 
obviously have some Coulomb disorder because of the 
inevitable presence of charged impurities. At low enough 
Ep (i.e. low enough doping), this Coulomb disorder is al¬ 
ways relevant for transport since it leads to zero conduc¬ 
tivity (i.e., infinite resistivity) whereas any short-range 
disorder (even if it is important at high carrier density 
for large Ep) becomes completely irrelevant at the Dirac 
point (as Ep —>■ 0) since it gives rise to finite conductivity. 
This implies that transport properties at the Dirac point 
{Ep = 0) will be completely determined by Coulomb dis¬ 
order leading to zero conductivity at the Dirac point (ex¬ 
cept for any finite-temperature effects we discuss above). 
This remains true, within the Boltzmann theory, inde¬ 
pendent of whether the DOS is finite or not at the Dirac 
point, and therefore, the Dirac point conductivity is al¬ 
ways likely to be semimetallic in nature (see Fig. 1) in 
the presence of any long-range disorder. 

A brief discussion of graphene (i.e., d = 2, a = 1) 
transport properties may be germane in this contexti*^ 
since d = 2 is the critical dimensionality for the 
semimetal-to-metal transition (for a = 1) implying that 
infinitesimal disorder produces finite DOS at the Dirac 
point in graphene, thus essentially making all undoped 
graphene samples everywhere being “metallic” with a 
universal Dirac point {Ep = 0) intrinsic conductivity 
of ap = Ae’^/irh. Such a universal conductivity has 
never been measured in graphene dc transport experi¬ 
ments where the transport measurements typically agree 
witb2ii2^ the Boltzmann theory predictions for long- 
range Coulomb disorder with cr ^ Ep oc n for higher 
carrier density with the conductivity crossing over at low 
carrier density {Ep ^ 0) to a broad nonuniversal mini¬ 
mum, which depends weakly on the amount of disorder 
(and other details), arising from inhomogeneous electron- 


hole puddle formation in the system which has been stud¬ 
ied extensivelyii2ri22iiMi2^ Our experience with graphene 
hints at what is likely to be seen in 3D Dirac transport 
when detailed experimental data become available. For 
Ep ^ kpT, we expect that the conductivity will fol¬ 
low the Boltzmann prediction^ going as Ep, and then 
depending on the details (e.g., temperature, amount of 
disorder, etc.), at some low values of Ep, there will be 
a crossover to the intrinsic a T'^ “semimetallic” be¬ 
havior which will persist as long as kpT > Ep. Eventu¬ 
ally, at some disorder-dependent low density, there will 
be a crossover to the puddle physics with the conductivity 
showing a nonuniversal broad minima around the Dirac 
point similar to what is extensively seen in graphene. 
This puddle physics is in some sense reminiscent of the 
physics of “rare regions” discussed in Ref. [l^ although 
the details are completely different since the puddles dis¬ 
cussed by us arise from the density inhomogenity inher¬ 
ent in Coulomb disorder in the absence of screening by 
carriers which must happen around Ep ^ 0^ The con¬ 
clusion, however, is the same that even the “semimetal¬ 
lic” phase at Ep = 0 always has a finite local density 
of states as the system nucleates electron and hole pud¬ 
dles (with fluctuating local density, but vanishing aver¬ 
age density for Ep = 0) in order to screen the impurity 
induced Coulomb potential. These finite local density 
of states at the Dirac point have directly been observed 
experimentally in graphene through imaging microscopy 
techniques and we believe that they must exist in 
3D Dirac systems also at low enough carrier doping^. 

One important caveat to note here is that all our trans¬ 
port results are obtained using the Boltzmann transport 
theory with the disorder-induced transport relaxation 
time calculated in the leading order Born approximation. 
Strictly speaking such a theory is, in principle, inappli¬ 
cable at the Dirac point with vanishing carrier density 
(where kp vanishes) since the dimensionless parameter 
kpl > 1, where I is the mean free path, must apply for the 
theory to be meaningful. In fact, this is a limitation not 
just for the Born approximation and the Boltzmann the¬ 
ory, but for all perturbative theoretical approaches to cal¬ 
culating the charge conductivity at the Dirac point since 
kpl is no longer a control parameter at the Dirac point. 
This is why our theory is carried out at finite carrier den¬ 
sity (i.e., finite kpl) and the Dirac point is approached as 
a limiting case of kp going to zero. As long as there is no 
singularity at the Dirac point, such a limiting procedure 
for obtaining the Dirac point conductivity should be at 
least qualitatively valid although our theory is strictly 
valid only for a low (but nonzero) carrier density so that 
kpl remains finite. We mention here that simple improve¬ 
ments of the theory, for example, using the self-consistent 
Born approximation or including the higher-order Born 
scattering terms is also not justified since all of these 
approximations are expansions in 1/kpl, which becomes 
ill-defined at the Dirac point since kp = 0 there! In par¬ 
ticular, we do not believe that the self-consistent Born 
approximation, which is sometimes used, is a better ap- 
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proximation at all since it is completely uncontrolled (and 
thus, its domain of validity unknown) and includes only 
selective higher-order impurity scattering diagrams leav¬ 
ing out many terms in each order in an ad hoc manner. 
One great advantage of the Born approximation-based 
Boltzmann theory is that it is known to be valid and accu¬ 
rate at higher carrier density (and has had great success 
in describing carrier transport properties in semiconduc¬ 
tors, metals, and graphene^), and thus, unless there is a 
quantum phase transition with decreasing carrier density, 
the theory should remain qualitatively valid as the Dirac 
point is approached by decreasing kp albeit with increas¬ 
ing quantitative inaccuracy. At this stage, we do not see 
a way to develop a better transport theory precisely at 
the Dirac point (i.e. kp = 0), which in any case is a 
set of measure zero not particularly relevant for experi¬ 
ments (e.g., any finite temperature immediately relaxes 
the zero-density condition, making the Boltzmann theory 
applicable in a formal sense). As we show in our work 
(see Fig. 1), the important physics is that the intrinsic 
(extrinsic) semimetallic (metallic) behavior manifests it¬ 
self in the charge conductivity for T > (<) Tp in gapless 
systems, and the strict zero-density or zero-temperature 
limit is unnecessary. 

We emphasize that the transport theory used in the 
current work is the Boltzmann theory where the effect 
of disorder is included only in calculating the transport 
scattering rate (and not in the DOS itself), and our 
goal is to correlate the Boltzmann transport (not crit¬ 
ical transport) behavior with the non-interacting DOS 
in the clean system. This approach is known to work 
very well in graphen o^^i^°i^^i^^i^^ , where the Boltzmann 
theory based calculated conductivit} d^i^°i^'^ agrees well 
with experimental results2ii2^i2^ for systems with arbi¬ 
trary disorder in spite of graphene being the critical case 
(satisfying d = 2a). If quantum critical properties be¬ 
come important in an experimental situation (and we 
know of no direct experimental evidence for the predicted 
quantum critical behavior dc = 2a in any physical sys¬ 
tem), then our Boltzmann theory is restricted to weakly 
disordered systems for d < 2a, and to systems with dis¬ 
order strength less than the critical value for d > 2a. 
But, our main goal of connecting the Boltzmann conduc¬ 
tivity with the clean system DOS for a gapless system is 
completely independent of these disorder-driven critical 
properties considerations. 

We conclude by saying that we have analytically ob¬ 
tained the scaling dependence of the Boltzmann con¬ 
ductivity a{Ep,T) in gapless systems with arbitrary 
band dispersion and dimensionality due to scattering by 
Coulomb and short-range disorder. Our results are gener- 
ically valid for both chiral and non-chiral gapless systems, 
and illuminate the behavior of semiclassical conductivity 
in the semimetallic (zero density of states) and the metal¬ 
lic (finite density of states) gapless systems (when the 
chemical potential is near the band touching or “Dirac” 
points), showing that the power law dependence (neglect¬ 
ing logarithmic corrections) of the disorder-limited con¬ 


ductivity could, in principle, distinguish the semimetallic 
behavior (cr ~ T'’') at high temperatures (T > Tp) from 
the metallic behavior (cr ~ Ep) at low temperatures. 
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Appendix A: Temperature dependent screening in 
the semimetal 

In deriving the main results of this paper (at finite tem¬ 
peratures), we neglected any temperature dependence 
arising from the screening itself, concentrating instead on 
the universal temperature dependence of the conductiv¬ 
ity arising from the energy averaging which dominates 
the high-temperature conductivity for T > Tp in the 
semimetallic phase (where the carrier density and hence 
the Fermi energy are small by definition). In principle, 
however, screening itself could have temperature depen¬ 
dence which can be quite subtle and rather important 
for non-chiral systems, particularly in 2D systems^^ with 
parabolic dispersion where the 2kp ‘kink’ in the zero- 
temperature polarizability leads to nonanalytic correc¬ 
tions of 0{T/Tp) in the Boltzmann conductivity aris¬ 
ing from screened Coulomb disorder—. In chiral systems 
of interest in the current work, 2kp back-scattering is 
suppressed, and therefore, the temperature dependence 
of conductivity is much weaker for screened Coulomb 
disorderi^i^ Nevertheless, the screening wave vector Qs 
may itself pick up a temperature dependence which we 
consider below (this effect is neglected in the main body 
of our work, but is straightforward to include^i^). 

At first, it may seem that since Qs is basically the 
density of states at the Fermi level, it cannot have any 
temperature dependence since the DOS, by definition, is 
temperature-independent. This is, however, untrue since 
at finite temperature, even for Ep = 0 at T = 0, an av¬ 
eraging over energy must be done in order to calculate 
any physical quantity. This averaging over the Fermi dis¬ 
tribution function leads to a T-dependence in screening 
even in the semimetallic phase at the Dirac point (and 
Qs is no longer just simply the DOS at Ep). 

The temperature dependent TF screening wave vector 
can be calculated from the d-dimensional polarizability 
at g = 0, 

n(9 = Q,T) = gJ dED^E) ■ (Al) 

Then, we have the Qs as 

g«(T)ocT^^. (A2) 
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For the linear dispersion (a = 1) we find that qs increases 
linearly with temperature regardless of dimensionality. 
We note that, for d > a, which is true for 2D and 3D 
Dirac systems, screening at the Dirac point is enhanced 
(i.e. qs increases) with increasing temperature, and thus 
leads to insulating or semimetallic temperature depen¬ 
dence (i.e. an increasing conductivity with increasing 


temperature) which is the same behavior as that from en¬ 
ergy averaging studied in the main part of the paper. For 
finite doping with Ep > kpT, the temperature depen¬ 
dence of qs is exponentially suppressed by exp(—T/Tp) 
and is no longer important. Also, disorder itself may sup¬ 
press the temperature-dependence of screening by mod¬ 
ifying the polarizability function as has been considered 
in the literature for parabolic systems22i2ii^. 
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